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The purpose of this note is to argue for the appropri- 
ate notion of stability in a theory with broken Lorentz 
symmetry that supports modes with phase velocities dif- 
ferent from the speed of light. In particular we are mo- 
tivated by the example of Einstein-aether theory, but 
our considerations are quite general. More specifically, 
we shall argue that the stability criteria imposed on this 
theory in Ref. [l[ are overly restrictive. The conclusion 
is that the theory is actually stable for an open set in the 
four dimensional coupling parameter space rather than 
for only a one dimensional subspace. The issue of sta- 
bility in Lorentz violating theories was also addressed in 
Refs. [1, Hj], which include arguments closely related to 
those advanced here. 

Einstein-aether theory is an example of a theory where 
Lorentz symmetry is dynamically broken. Aside from 
matter, the fundamental fields are the spacetime metric 
g a b and a timelike unit vector field u a , the "aether" (see 
Ref. [H for a review). Flat spacetime with a constant 
aether is a solution to the theory, and linearized pertur- 
bations of this solution satisfy second order hyperbolic 
equations. There are modes with five different polariza- 
tions: two spin-2, two spin-1, and a single spin-0 mode. 
For all these modes, the frequency u> and spatial wave 
vector k defined relative to the rest frame of the aether 
satisfy a gapless dispersion relation, u> 2 — vfk 2 , where i 
labels the spin. The squared velocities vf depend on the 
coupling parameters in the Lagrangian, and are generally 
different from each other and different from the "speed 
of light" c defined by the null cone of the metric g a b. The 
conditions v 2 > impose inequalities on the coupling pa- 
rameters, guaranteeing that the frequency is real, so the 
perturbations do not grow exponentially in time if the 
spatial wave vector is real [5]. Another set of inequal- 
ities implies that the energy carried by these modes is 
positive 

urn, 

and these inequalities can be satisfied si- 
multaneously with the former stability inequalities. (As 
of yet, no nonlinear extension of this positive energy re- 
sult is known, except in the special case of static spherical 
symmetry Q.) 

It was recently argued in Ref. [l[ that one should re- 
quire the above stability criteria, i.e. real frequency and 



positive energy, not only for modes with a real wave vec- 
tor in the aether frame, but more generally for modes 
with real wave vector in any Lorentz frame defined with 
respect to the metric g a b, and for energy defined in any 
Lorentz frame. In a non-Lorentz invariant theory this 
is obviously a much stronger requirement, and in fact it 
was concluded that Einstein-aether theories are unsta- 
ble except for a small number of special cases in which 
all modes propagate at exactly the speed of light 1 . Wc 
shall now argue, however, that these stronger conditions 
are not required by stability of the theory, and are not 
justified given the structure of the theory. 

In fact the reasoning of Ref. [l[ applies to any linear 
theory with modes propagating at different speeds, not 
only to Einstein-aether theory. Also, the dynamics of the 
metric and aether themselves play no essential role except 
to define those modes. Hence we will discuss the simpler 
setting of fields on a spacetime with a fixed Minkowski 

metric rj a b (with signature (H )) and a fixed timelike 

unit vector u a . A free scalar field ip that propagates at 
speed v with respect to the rest frame of u a is minimally 
coupled to the effective (inverse) metric 

S^A' + ^^-uV), (1) 

with Lagrangian density 

C = | V=*7 d W = | (d?<p - v 2 df^j , (2) 

where the second expression is written in the Minkowski 
coordinate system (t, x 1 ) of the metric r) a j,, adapted to the 
rest frame of u a . We consider this model with arbitrary 
positive values of v. 

In the context of Einstein-aether theory, there is good 
reason to allow v to be greater than c. If the coupling con- 
stants are chosen so that the post-Newtonian preferred 
frame parameters of the theory are in agreement with 
observational constraints, then the positivity of the en- 
ergy (in the aether frame) requires v > c. Also, to satisfy 
the vacuum Cherenkov constraint for ultra high-energy 
cosmic rays, all v > c are allowed, but any v less than c 
must be extremely close to c [1]. 



1 Actually, in Ref. []J the decoupling limit was taken. That is, the 
'Electronic address: wdonnell@umd.edu metric is fixed to the Minkowski metric and not varied in finding 
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The dispersion relation for a scalar field with La- 
grangian density © is w 2 = v 2 k 2 when expressed in 
terms of components of the wave 4-covector k a in the 
aether frame. More precisely, the wave phase is k^x^ — 
uit + kiX 1 , and k 2 = kiki. For real spatial wave vectors 
hi, the stability requirement that the frequency be real 
amounts to the condition v 2 > 0. This condition guar- 
antees that any solution that is a superposition of plane 
waves on a constant t surface is stable. 

One of the further stability criteria of Ref. [l| is the 
demand that the frequency be real also for plane waves 
on any Lorentz-boosted constant time surface. To deter- 
mine what that implies, we may reexpress the dispersion 
relation in terms of the components of the wave 4- vector 
in the boosted frame as follows. 

The metric r\ a b can be used to define a set of frames, re- 
lated in the usual way by Lorentz transformations. With 
respect to such a frame moving with velocity /3, the new 
time and space coordinates are given by 

t' = -y(t-0x\\), (3) 
x\ = l{x\\-fc), (4) 
x'x = xx, (5) 

where || and j_ refer to the components parallel and per- 
pendicular to the boost direction, and we use units with 
the metric speed of light equal to unity, c = 1. The frame 
velocity will be taken to be positive, and is assumed to 
be less than the speed of light, < (3 < 1. 

The covariant (as opposed to contravariant) frequency 
and wave 4-vector components in the boosted frame are 
given by 

J = 7 (w + /3fc||), (6) 
fcfl = 7(fc||+0w), (7) 
k' x = fc x . (8) 

The dispersion relation in terms of these boosted compo- 
nents takes the form 

(1 - v 2 /3 2 )uj' 2 + 2/3k\ { (1 - v 2 )uj' 

+ {P 2 -v 2 )k l2 -v 2 (l-p 2 )k' 1 2 = 0, 

where we have multiplied by a factor (1 — (3 2 ) = 7~ 2 for 
convenience. This is a quadratic equation for a/, so the 
roots are real for real k' if and only if the discriminant is 
positive, 

v 2 (l - f3 2 ) 2 k\ { 2 + v 2 {l - v 2 p 2 ){\ - /3 2 )fc ± 2 > 0. (10) 

Since f3 < 1, this can be negative only if the term (1 — 
v 2 f} 2 ) is negative, which occurs only if v > 1 and (3 > 1/v. 
We note that this is just the condition for the constant 
t' surfaces to be timelike with respect to g( v ) a b ©■ The 
frequency uj' then has a nonzero imaginary part when 
k' ± /k'^ is sufficiently large. 

Thus for v > 1 and /3 > 1/v there exist solutions with 
real wave vectors and complex frequencies in the boosted 



frame. Such modes grow exponentially in the time coor- 
dinate t' of that frame. Whether or not this indicates an 
instability comes down to the question whether or not 
these solutions are part of the physical phase space of 
the theory. 

As pointed out in Ref. the wave vector kn — 
7(fef, —j3u>') in the rest frame of the aether will be complex 
for such modes, so on a constant t surface the solution will 
blow up exponentially at spatial infinity. These solutions 
therefore do not satisfy the usual boundary conditions 
that define the phase space of the theory on the constant 
t slices. A consistent theory can be defined by adopting 
a regular boundary condition on the constant t slices, ex- 
cluding these solutions. This is what is ordinarily done 
in a Lorentz invariant theory. For example, one could 
require that the solutions have compact support, or that 
they be Fourier transformable on those slices. Moreover, 
since these boundary conditions are preserved by t evolu- 
tion, the theory so defined preserves the time translation 
symmetry of the background. Also, one would define the 
same phase space imposing these boundary conditions on 
any other surface that is spacelike with respect to g^ v ) a b- 

One may ask whether a consistent theory could instead 
be defined by adopting a regular boundary condition on 
the constant t' slices. If so, this would raise the ques- 
tion of which is the correct phase space. But it appears 
that this can not be done in a natural way. For certain 
regular initial data on a given constant t' surface the 
corresponding solution grows exponentially with t' . In 
any other frame this solution will contain complex wave 
vectors and will therefore diverge asymptotically on the 
constant time slices of that frame. This means that, un- 
like the case for surfaces which are spacelike with respect 
to 

9{v)ab, the phase space defined by regular data on con- 
stant t' surfaces is different for every value of f3 greater 
than 1/v. 

Moreover, for /? greater than 1/v, the phase space ob- 
tained by requiring regular initial data on a fixed t' slice 
depends not only on the particular value of /3, but on the 
particular choice of slice. For example, suppose the ini- 
tial data on a particular surface t' = t' possesses a well- 
defined Fourier transform. In Fourier space, the wave 
equation then reduces to an infinite number of uncou- 
pled ordinary differential equations that may be solved 
to obtain the Fourier transform of the solution on a dif- 
ferent slice t 1 = t[. For modes with sufficiently large 
kj_, the solutions to the differential equation grow expo- 
nentially with k±. The solution at t[ in Fourier space 
therefore does not in general have a convergent inverse 
Fourier transform. Such initial data on the t' surface 
do not correspond to any choice of initial data on the t[ 
surface. This means that the phase space defined in this 
way depends not only on the choice of time coordinate 
t' but also on the arbitrary value t' of that coordinate, 
breaking the time translation symmetry of the theory. 

Another reason to reject a "t'-phase space" formulation 
is that allowing for arbitrary initial data at t' = t' is 
unjustified in the context of a causal theory in which 
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the ip field interacts with other degrees of freedom. A 
simple way to see the problem is to allow for an external 
source term in the field equation for ip. One can then ask 
whether the source could generate data at t' Q that would 
lead to an exponentially growing solution. The answer 
is no unless (perhaps) if the source is turned on in the 
infinite past. As explained above, any such solution will 
blow up exponentially at spatial infinity on all constant 
t surfaces. If the source is turned on at a finite time, its 
effects cannot propagate any faster than v in the aether 
frame, and so the solution can not blow up at spatial 
infinity at any finite time. 

The preceding argument depends on a choice of bound- 
ary condition for the solution generated by the source, 
which is equivalent to the choice of Green's function for 
the wave equation. We implicitly adopted the retarded 
Green's function, which vanishes for t < 0. One might 
ask whether the argument would continue to hold us- 
ing a i'-retarded Green's function that would vanish for 
t' < 0. It appears, however, that no such Green's func- 
tion exists for f3 > 1/v. This can be seen as follows. 
A standard method for constructing Green's functions is 
via the Fourier transform 

/ p ik' -x'~- iuj' t' 
<Pk'dJ — -— (11) 
(u)' - - 

where oj'± are the roots of the dispersion relation ([9]). 
This integral can be performed along any contour that 
begins at — oo and ends at +oo along the real axis. If 
the integral converges, then the wave operator acting on 
G(t',x') can be moved inside the integral, canceling the 
denominator. The remaining integrand has no poles, so 
the contour can be freely deformed to lie along the real 
axis, yielding a representation of the Dirac delta function. 
To obtain the retarded Green's function, the u>' integral 
is performed along a contour that passes above all the 
poles, so that for t' < the integral vanishes. For t' > 0, 
both poles are enclosed by the contour that can be closed 
in the lower half plane, and the u> integral yields 

d 3 k , e ik .» e e 

- u'_ 

For large k' ± , the roots behave as u)'± ~ ±ik' ± , so that the 
integrand in (|12j) grows exponentially with k'^ and the 
integral does not converge. A i'-retarded Green's func- 
tion therefore cannot be found by this standard method, 
which strongly suggests that such a Green's function does 
not exist 2 . 

A second sign of possible instability discussed in 
Ref. [l[ is that the Hamiltonians generating t' trans- 
lations can be unbounded below whenever v =/= 1. In 



2 Note that the exponential instability alone does not account for 
the absence of a retarded Green's function. For example, in the 
case of a tachyonic scalar field with a negative m 2 , the instability 
occurs only at low k, so the convergence of the Green's function 
is not spoiled, and a retarded Green's function exists. 



particular the Hamiltonian of linear perturbations is un- 
bounded below precisely when /3 > v ot /3 > 1/v (the con- 
dition in Ref. [l| was expressed in terms of the coupling 
constants of the theory rather than the mode speeds, nev- 
ertheless the two conditions are equivalent). This can be 
understood in a simple way as follows. 

Let £ a denote the t' translation 4-vector. The Hamil- 
tonian generating £ a translations can be written as an 
integral over an initial data surface E, 

H 6 = [ T b a en b d 3 Z, (13) 

where T b a is the canonical energy-momentum tensor 

BC 

^ Tb ^Ww) da ^ csl (14) 

rib is the unit normal covector, and <i 3 E is the surface 
volume element, both normalized with respect to 7y ao . 
Positivity of is ensured when T h a £ a nfc is positive. In 
the case v < 1 the ip field is "subluminal" relative to 
rjab, so for j3 > v the vector £ a can be timelike relative 
to rjab but spacelike relative to the effective metric g( v ) a b 
for the ip field. In this case H$ is in effect a component 
of the momentum, not the energy of tp, which is clearly 
not bounded below. In the case v > 1, for /? > 1/v the 
t' translation vector £ a remains timelike with respect to 
9(v) ah j but the constant t' surface becomes timelike with 
respect to g( v ) a b- In this case is the flux of energy 
through a timelike surface, and is no longer expected to 
be bounded below. Moreover, if the surface E is timelike 
with respect to gt v ) a b there is no guarantee that is 
conserved under t' translation, because the current T b a £ a 
can flow out through the boundaries. 

In conclusion, while we take no issue with the computa- 
tions of Ref. [l[ , the inference of instabilities in Einstein- 
aether theory when the mode velocities differ from c is 
unwarranted. A proper identification of the phase space 
of the theory eliminates the exponentially growing so- 
lutions. The Hamiltonians that were found to be un- 
bounded below actually correspond either to momenta 
or to energy fluxes across timelike surfaces. The Hamil- 
tonian generating time translations in the aether frame is 
bounded below and plays the usual role of the energy in 
governing stability. It is therefore sufficient for stability 
to impose the conditions of real frequencies and positive 
energy in the aether frame. The opposite conclusion was 
reached in Ref. [l| by considering the Lorentz symmetry 
of the background metric to be a physical symmetry of 
the phase space of linear perturbations. Since the back- 
ground aether breaks this symmetry, that viewpoint is 
untenable. 
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